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Making use of the Bogoliubov-de Gennes equation, we study the quasi-particle spectrum 
and the vortex core structure of a single vortex in quasi 2D s-wave superconductors for 
small PfCo, where pf is the Fermi momentum and = vf/Aq is the coherence length 
{h = 1). During our numerical calculation the particle number is conserved for each pf^o- 
In particular, we find that there are only 1 or 2 bound states for pf^o ~ 1- Also, for pfCo = 1> 
the Kramer-Pesch effect ceases to exist at around T/Tc — 0.3. 

§1. Introduction 

There are renewed interests in vortex structure since the discovery of high T(. 
cuprate superconductors. It is possible that the sunecconductivity of high Tc cuprates 
can be chaxacterized as d-wave superconductivity H^'cP and is very close to the quan- 
tum limit. B^'Q^ 

Schopohl and Maki previously studied the quasi-particle spectrum around a 
single vortex line in terms of a quasi-classical equationQ^'Q) and predicted a clear 
four-fold symmetry for a d-wave superconductor. Ix Then, a beautiful STM study 
of vortices in YBCO monocrystals was reported. There are several interesting 
results of these studies. First, it was found that the quasi-particle spectrum exhibits 
circular symmetry; there is no trace of four-fold symmetry Second, there appears 
to be only a single bound state with energy ^ jA{0), where A{0) (= 260 K) is 
the superconducting order parameter at T = K. A similar bound state energy 
was observed earlier by a far-infrared magneto-transmission from a YBCO film by 
Karrai' et al.EF According to the analysis of bound states around the vortex core by 
Caroli and co-workers, E^P there should be thousands of bound states. Of course, in 
usual s-wave superconductors we have pf(,o ~ 10'^ ~ 10^, where pp is the Fermi 
momentum and = vp/^iO) is the coherenca length (h = 1). Thus perhaps the 
single bound state in YBCO suggests pfS,o ~ 1.3* At first sight this suggestion seems 
outrageous, since this implies Ep ~ 200 ~ 500 K in YBCO. By analyzing the spin 
gap at T = K observed in YBCO monocrystals with inelastic neutron scattering 
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by Rossat-Mignod et al.,lllP we can deduce chemical potential /i given asO 

= -345(x - 0.45) K, (M) 

where x is the oxygen dopage corresponding to YBa2Cu306+x- Then, for optimally 
doped YBCO we obtain ^ = -190 K. 

In Dottinger et al.,u^ the flux-flow resistance of 60 K YBCO measured by Mat- 
suda et al.EJ) is analyzed, and they identifiecLthe Kramer- Pesch effect characteristic 
of a superconductor in the clean limit. EJ)'tM) On the other hand, apparently the 
Kramer-Pesch effect is absent in 90 K YBC0,ll3) implying again thar perhaps 90 
K YBCO is in the quantum limit. Here, the Kramer-Pesch effect causes the vor- 
tex core to shrink with decreasing temperature due to the decrease in the occupied 
bound states around the vortex line. The core size can be expressed as 

This reduction of the core size results in nonlinear conductivity. Ichioka et 

al.c3 found the Kramer-Pesch effect in d-wave superconductivity with the help of a 
semi-classical approach developed in Ref. 0). 

But if pf(,o ~ 1) this implies that the semi-classical approach introduced in 
Refs. ^) and^) is no longer reliable in studying the vortex in high cuprates. For this 
reason, Morita, Kohmoto and Makit3* studied the Bogoliubov-de Gennes equation 
for a single vortex in d-wave superconductors. Indeed, choosing pf(,o = 1-33, they 
were able to describe gross features of the STM result for YBCO. On the other hand, 
they have not attempted the self-consistent calculation due to numerical difficulties. 
Nevertheless, they discovered that there is a single bound state for pfS,o = 1-33 in 
the vortex of d-wave superconductors. Further, there are low energy (E < 0.1Z\) 
extended states with four legs stretched in the four diagonal directions (±1^1,0). 
Recently, these extended states were rediscovered by Franz and Tesanovict^ in a 
model somewhat different from that in Ref. |l^. Further, Franz and Tesanovic 
claimed there should be no bound states, based on a different model which contains 
at least 3 arbitrary parameters. We believe that the absence of bound states is 
either due to a) neglect of particle number conservation or b) the rather strcmg 
Coulomb repulsion they introduced. It is rather surprising that Yasui and Kitali^P 
and Takigawa et al.S obtained similar results as Ref. [l8| ). Clearly in these later 
works, particle number conservation has been neglected. On the other hand, we have 
shown recently there are a few bound states within the weak coupling model, as in 
Ref. ^), down to Pf£,o — 1.^^ Also, this shows clearly that number conservation 
is crucial, as pointed out by van der Marel.S Also, the above results confirm the 
validity of the work in Ref. ^) rather than that in Ref. 18). 



In this paper we study the quasi-particle spectrum around a single vortex in 
s-wave superconductors in the quantum limit. This problem was previously consid- 
ered by Hayashi et al.c2P However, as we have mentioned, they ignored the number 
conservation, which is very crucial. For example van der Marel has shown that the 
chemical potential depends on the temperature if number conservation is imposed. 
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Further, in Ref. 23), the energy cutoff is made arbitrarily, so that their result is not 
reliable, as we will show. Also, as we will show, the r-dependence of A(r) near 
r = is completely determined by UQ{r)vo{r), where uo{r) and VQ^r) are spinor wave 
functions for the lowest bound state. 

In the following we concentrate on the cases pf^o = 1)2 and 4 for simplicity, 
and study the quasi-particle spectrum and the shape of |^(r)| as a function of 
temperature. A preliminary result on this has appeared in Ref. p^). 

§2. Bogoliubov-de Gennes equation 

The spatial dependence of order parameters in superconductivity is described 
by the Bogoliubov-de Gennes equation. We consider the case at a small magnetic 
field near Hd, so that there is a single vortex. Therefore we can ignore the vector 
potential. In this case, the Bogoliubov-de Gennes equation becomes 

^ - fi] Un (r) + A (r) Vn (r) = E^Un (r) , (2- la) 



2m, 
1 



2mp 



V2 - /i U„ (r) + A* (r) un (r) = E^Vn (r) , (2- lb) 



where Un {t) and Vn {r) are quasi-particle wave functions. 

We take the z-axis to be along the vortex line. We consider the nearly two- 
dimensional case for simplicity, where the kinetic term associated with the z direction 
is negligible. In the following, we merely consider the two-dimensional case, and we 
use cylindrical coordinates. Taking the gauge as A[r) = |Z\(r)|e~!^ the angular 
momentum of each eigenstate becomes half of an odd integer, m + ^.113 Then Unir) 
and Vn{r) become as follows; 

Unir, 0) = Unm{r)—=, (2-2a) 
v 2tt 



gi(m+l)e 

Vn{r,9) = Vnmir) — 7=^- (2-2b) 
V ZTT 

Following Gygi and Schliiter,ii' we apply the Fourier-Bessel expansion with basis 

\/2 T 

(t>mj{r) = — zJm{ajm-^) (2-3) 

to 

Unmif) arid Vnmif): where Oijm is the j-th positive zero of the Bessel function of 
m-th order, Jm{x). Then the wave functions become as 



Unmir) = ^Unmj4'mj{r), (2-4a) 
j 

Vnmir) = Vnmj4>m+ljir). (2-4b) 



J 



Here, the boundary condition is such that wave functions are zero at the edge of the 



4 



M. Kato and K. Maki 



disk with radius R. Then the Bogohubov-de Gennes equation becomes 



1 



2mp 



R 



1 



2me 



a 



j m+l 

R 



imj + ^ 



nmj ) 



(2-5a) 
(2-5b) 



where /\jij2 is given as 
Z\ 



<Pmh{r) \^{r)\ (pm+ij2 {r)rdr. 



(2-6) 
(2-7) 



The order parameter is given as 
\^{r)\=9 E EE' 

I ^nm 

where g is the interaction constant, Ec is the cutoff energy, and f{E) is the Fermi 
distribution function. 

Also we impose particle number conservation. 



N 



2 lY.^^^'nir)? f{En) + \vn?{l- f{En))}d^r, (2-8a) 

•' n 

2EE[k 



f{En 



[l-fiEn)) 



(2-8b) 



where is the total number of particles. The chemical potential is determined by 
this equation. Here is the particle number of the normal state with specified pp 
at T = 0. We solve these equations, Eqs. (2-5), (|2-7D and (2-8), self-consistently. We 
fix the cutoff energy (Ec = 5Z\o for the case of largest Fermi momentum) and the 
interaction constant, and choose the Fermi momentum pp so that pf£,o = 1)2 and 4. 

§3. Results 

In the numerical calculation, the radius of boundary is taken as ii = lO^O) and 
maximum angular momentum and number of zero points of the Bessel functions are 
taken so that all of the quasi-particle states within the cutoff energy are taken into 
the calculation. 

Before proceeding to the numerical results, here we discuss the temperature 
dependence of the chemical potential for the superconducting state. The chemical 
potentials for the normal state (jJm) and the superconducting state (/x^) of the two- 
dimensional system are given asc3 



pfin (T) = ln (e^A*"(o)-l) 



2/x„ (0) = lis (T) +E,- sjE^ + A {Tf + {TY + A (T) 
- In + - In h + e"'^^^ 



l + e 



(3-1) 



(3-2) 
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The temperature dependences of the chemical potential for pp^^ = 1 and 4 are shown 
in Fig. 1. As pointed out by van der Marel,E2P the chemical potential decreases with 
decreasing temperature in the superconducting state. This is opposite to the behav- 
ior in the normal state. As we can see from the figure for pp^o = 1; the chemical 
potential becomes very small in the superconducting state, and the chemical poten- 
tial of the normal state drops significantly from T = to T = Tc. For larger pp, this 
temperature dependence of the chemical potential becomes much weaker. There- 
fore the temperature dependence and the pp dependence of physical quantities may 
change from those with the constant chemical potential calculation in the quantum 
limit. 




Fig. 1. Temperature dependence of the chemical potential for pf^o = 1 (a) and pf£,o ~ 4 (b). The 
solid curve is the chemical potential for the superconducting state, and the broken curve is that 
for the normal state. 



3.1. Temperature dependence of order parameter 

The temperature dependence of the order parameter is shown in Fig. 2. Near 
the boundary r/^o = 10, the order parameter goes to zero, and there is a peak before 
that. These are effects from the boundary condition and the finite size. But the core 
structures are not affected by the boundary condition. 

For pp^o = 4 at low temperature {T/Tc < 0.1) there is a shoulder in 
at the vortex core. This feature comes from bound states in the vortex core. This 
can be seen from Fig. 3, where contributions of the scattering states and the bound 
states are shown separately. 

For pp(,o = 1, the peak position of the contribution to the order parameter 
is located slightly outside of the vortex core. Therefore the core structure is not 
so strongly affected by the bound states. For larger pfCo(~ 16), Hayashi et al.B5' 
demonstrated the oscillation of the order parameter on the inside of the vortex core. 
They argued that this oscillation is Friedel oscillation, but from our figures it appears 
that the origin of the oscillation is the discreteness of the bound states. 

Slightly increasing the temperature, the bound state contribution decreases 
rapidly, and the scattering state contribution remains almost the same. In this 
temperature range, the core size is dominated by the bound states. Above this 
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Fig. 2. The temperature dependence of the or- 
der parameter A{r) for (a) pf^o ~ 4, (b) 
Pf^o = 2 and (c) Pf^o = 1- The order pa- 
rameter is normalized with Aq = A{r = 
5.5^0,7" = 0) for each pfCo- 
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Fig. 3. Contributions from the bound states 
and the scattering states to the order pa- 
rameter for (a) pfCo = 4, (b) pf^o = 2 
and (c) PfCo ~ 1, where Aq — A{r = 
5.5Co,r = 0). 



temperature region, the scattering state contribution decreases with increasing tem- 
perature, which is the same behavior as that of the order parameter of the uniform 
solution. 

We also plot the quasi-particle wave functions Un^{r) and Vn^{r) of three bound 
states ior pfS,o = 4 in Fig. 4. From these figures, it can be seen that main contribution 
to the core structure comes from the lowest energy bound state. 

Also, we can see that u{r) of the lowest energy bound state behaves like an 
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s-wave function, and v{r) behaves like p-wave function. Similar behavior also can 
be seen for the second and third lowest bound states. They belong to the angular 




Fig. 4. The quasi-particle wave functions u{r) 
and v{r) and their product with a thermal 
factor for (a) the first, (b) the second and 
(c) the third of the lowest energy bound 
states. 
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Fig. 5. The quasi-particle spectrum of the vor- 
tex state (O) and the uniform state (dashed 
curves) for (a) Pf^o = 4, (b) PfCo = 2 and 
(c) pfCo = 1. 



momentum m and m -|- 1 states, respectively. 
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(b) 1 



3.2. Quasi-particle spectrum 

In Fig. 5 we give quasi-particle energies with angular momentum m + ^, and 
we also plot those of the uniform state in order to compare the vortex state and the 
uniform state. For small angular momentum, the lowest energy states become bound 
states, lowering their energies and slightly increasing the energies of higher energy 
states. For larger angular momentum, the quasi-particle energy is almost same for 
vortex and uniform states, and there is no bound state. 

Although it is difficult to determine the 
boundary of these two cases in the an- 
gular momentum, the number of bound 
states may be counted as 18 for pfS,o = 
4, 5 for pf(,o = 2, and 1 or 2 for pfCo = 
1. Compared with our previous calcula- 
tion without number conservation, the 
number of the bound states decreases. 
This effect is clearer for smaller pf^o- 

Such a relation between the num- 
ber of the bound states and Fermi mo- 
mentum is important for d-wave super- 
conductors, because the local density of 
states changes greatly due to this rela- 
tion. Therefore, it is worthwhile to point 
out this relation for the s-wave case, 
and in order to obtain the correct rela- 
tion, our number conserved calculation 
is needed. 

The temperature dependence of the 
quasi-particle energy is shown in Fig. 6. 
In this figure we plot the quasi-particle 
energies of several bound states normal- 
ized with the order parameter A{T) at 
r = 5.5^0- From this figure we can see 
that for smaller pf^o and higher energy 
states, the quasi-particle energy varies 
in proportion to the order parameter 
with temperature. But for lower en- 
ergy states with larger pf^o-, the quasi- 
particle energy increases rapidly with 
decreasing temperature and becomes 
constant at low temperature. This be- 
havior of the energy of the lowest bound 
state waajpreviously noted by Gygi and 
Schliiterc3 for much larger pf^o, but 
their result shows that the energy increases with decreasing temperature even at 




Fig. 6. The temperature dependence of the 
quasi-particle energies of the bound states 
in the vortex state. The energies are nor- 
mahzed by the order parameter A{T) at 
r/io ~ 5.5. (a), (b) and (c) are for pf^o = 
4, Pf€,o = 2 and Pf£,o ~ 1, respectively. 
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low temperature because of the large value of pf^o- 
3.3. Core radius 

As mentioned in 3.1. the core structure depends on the bound states. As defined 
by Kramer and Pesch,EJ) we calculate the core radius as 



1 

Ti 



A{r) 



lim ^ , 
r^o rZ\o 

g 



X! X! UnOjiVnOji [1 - 2/ {Eno)] 4>0ji (0) 
E„o<Ecjij2 



dct>ij2 
dr 



(3-3a) 
(0). (3-3b) 



We plot total and the contribution from the scattering states to (^^^ in Fig. 7. 

The temperature dependence of 
is almost linear in the intermediate tem- 
perature region, according to Kramer 
and Pesch.llS For pp^o = 4, there is a 
substantial decrease of the core radius 
at low temperature because of empty 
bound states. Although for smaller pp^o 
there is a bound state contribution, the 
peak position for its contribution to the 
order parameter is nearly at the edge 
of the core. Therefore the core radius 
is not strongly affected. But for large 
Pf^o, the peak of the contribution of the 
bound states is well inside the core, so 
the core radius shrinks significantly. 

Our result exhibits saturations of 

at T/Tc w 0.1 for pp^Q = 4, at 
T/Tc « 0.2 for pFio = 2 and at T/T^ ? 





p^=4 total 




. p£„=4 scattering state - — 




\ Pp^°=2 total 




\ p £ =2 scattering state - 




\ p^^°=1 total 




\^Pi-^j,=1 scattering state " 



















0.2 



0.4 

T/Tc 



0.6 



0.8 



Fig. 7. The temperature dependence of 

for pf£,o = 1,2 and 4. Also contributions 
from scattering states for each pfCo are 
plotted. 

0.3 for pf^q = 1- This saturaion comes 



from the disappearance of the occupied bound states at low temperature. In Fig. 3 
in Ref. 23), .^i is found to saturate at T/Tc « 0.1 for small pf^o- However, the au- 
thors of that work did not impose the particle number conservation condition, and 
therefore the temperature dependence of .^i is different from that in our result. 

3.4. Current density and magnetic field 
The current density is calculated from 



2mpi 



{/(^nm)<^(r-)Vn„„(r) + [1 - r;„„,(r)V<^(r) - h.c.} . 

nm 

(3-4) 



There exists rotational symmetry, so the current has only a 6 component: 

je{r) = — \f{Enm)— \Unm{r)\^ - [1 - f{Enm)] \Vnm{r)\^ 



(3-5) 
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(a) (b) (c) 



1.2 




01 234501 234501 2345 

r/^o r/^o r/^o 

Fig. 8. The current density at T = O.lTc for (a) pp^o = 4, (b) pF^o = 2 and (c) pp^o = 1- The 
current density is normahzed by its maximum value. The order parameter is also shown for 
comparison. 




Fig. 9. The spatial dependence of the mag- 
netic field H a.tT = O.lTc for pp^o = 1,2 
and 4. 



We show this current density at 
T = O.lTc in Fig. 8 for each pf(,o- 
We see that the peak position of the 
current density is almost the same as 
the peak position of the lowest energy 
bound state. Therefore, for a smaller 
value of pf(,o, the peak is located farther 
from the vortex core. Also, for smaller 
Pf£,0: the contribution of the scattering 
states is the same as that of the bound 
states, although for larger pfCo the con- 
tribution of the scattering states is op- 
posite to that of the bound states. Therefore there is no simple relation between the 
core size and the peak position of the current; that is, they are not exactly the same. 
From Maxwell's equations, we calculate the magnetic field parallel to the 2;-axis 

as 

H,{r) = — / je{r)dr. (3-6) 

C Jr 

We have normalized this with (/)o/27rA^, and show in Fig. 9, where A = {me(? / Anne^)^^'^ 
is the penetration depth and (j)Q = hc/2e is the flux quantum. From this figure, we 
find that for smaller PfCo, the distribution of the magnetic field is extended farther 
to the outside of the vortex core. 

3.5. Local density of states 

In order to compare with STM experiments, we calculate the local density of 
states with thermal average, 

N{r, E) = - Y.{\Unm{r)\''f{Enm - E) + \Vnm{r)\'' fiEnm + E))} . (3-7) 



We also take the spatial average with a Gaussian distribution with standard de- 
viation 0.1^0- This is shown in Fig. 10. In this figure there is oscillatory behavior for 
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IS 



the 



average 



E/Ao > 1.0, where 

value of the order parameter outside the 
vortex core. This feature results from 
the finite size of our system, and it is an 
artifact. In addition to this, there are 
several peaks due to the discrete bound 
states for energies less than the energy 
gap {E < Aq), and there is a particle- 
hole asymmetry, in contrast to the re- 
sult of Gygi_and Schliiter. i) Wang and 
MacDonaldE£P demonstrated this asym- 
metry by solviiLK the lattice model, ajid 
Morita et al.lilP and Hayasi et al.Ba) 
pointed it out previously. At pf(,o = 1 
there is only a peak near the energy gap, 
in contrast to the d-wave case, where for 
smaller pp^o there are only shoulders in- 
side of the energy gap.t 



§4. Conclusion 

In summary, we have solved the 
Bogoliubov-de Gennes equation with 
the number conservation condition for 
a single vortex state in s-wave super- 
conductors for small prCo, and we have 
obtained the electronic structure around 
the vortex core. The number con- 
servation condition strongly affects the 
quasi-particle spectrum and the core 
structure, because the chemical poten- 
tial varies strongly with the temperature 
in the quantum limit (pf'^o = 0(1)), Fig 
especially for smaller pf£,o- Also, we 
have elucidated the effect of the discrete 
bound states on the order parameter 
structure, the local density of states, the 
current density and the magnetic field. 
For smaller pf(,o, the contribution of the 

bound states to A (r) decreases, and the actual shape of A (r) depends on the peak 
positions of the product of the wave functions of bound states. Also, the discreteness 




f /A 

10. The local density ol states N{r, E) at 
T/T, = 0.1 for (a) ppio = 4, (b) pf£,o = 2 
and (c) pf£,o = 1- N{r,E) is normalized 
by the density of states of the normal state 
iVo. 



of the bound states is mori 
Recently, Sonier et al.l 



pparent for smaller pf^o- 

claimed that they obtained the core size of the vortex 
in NbSe2 by measuring the peak position of the current density. However, there are 
a few problems in their procedure. First, as we have shown, there is apparently no 
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simple relation between the peak in the current density and the core size. Second, 
even if such a relation did exist the relation between these quantities in the clean limit 
would be different from that in the dirty limit, as given by the result in Ref. 27). In 
particular, the core size in the clean limit is much smaller than that in the dirty limit 
with respect to the peak position in the current density. At least those authors should 
have analyzed their data in light of the clean limit calculation by Gygi and Schliiter. 
cB Third, the muon spin rotation experiment measures the average of the spatial 
distribution of the magnetic field. Therefore it is not suited to extract the current 
distribution. Rather, it should be used to extract the magnetic-field dependent 
magnetic penetration depth, which will provide more insight mto the vortex state. 

As compared with the result for d-wave super conductor s,clP our result exhibits 
systematic decrease in the number of bound states with decreasing Fermi momentum. 
In the d-wave case, there are bound states and extended states around a vortex. 
Therefore detailed study is needed for d-wave superconductors. Such a calculation 
is now in progress. 
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